Quasi-deterministic transport of Brownian particles in an oscillating periodic potential 
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We consider overdamped Brownian dynamics in a periodic potential with temporally oscillating 
amplitude. We analyze the transport which shows efTective diffusion enhanced by the oscillations 
and derive approximate expressions for the diffusion coefficient. Furthermore we analyze the effect 
of the oscillating potential on the transport if additionally a constant force is applied. We show 
the existence of synchronization regimes at which the deterministic dynamics is in resonance with 
the potential oscillations giving rise to transport with extremely low dispersion. We distinguish 
slow and fast oscillatory driving and give analytical expressions for the mean velocity and effective 
diffusion. 

PACS numbers: 



In recent decades there has been intense research on 
transport phenomena in Brownian dynamics. Numerous 
pubUcations analyze the diffusive transport in spatially 
periodic potentials and report interesting behavior of the 
diffusion coefficient [l|43| . Diverse sorts of ratchets have 
been investigated [^-|6[ where particles move in periodic 
potentials and time dependent forces or modulations of 
the potential enhance a directed flow. Most of the re- 
search focused on mean drift within such systems and 
on the conditions under which the directed transport is 
maximized. Also the diffusion coefficient in these tempo- 
rally changing potentials have been calculated [18[ and it 
was proposed to take the diffusion coefficient to evaluate 
the precision of stochastic directed transport [7|. 

In this spirit we will focus this work on overdamped 
Brownian transport in a temporally oscillating and spa- 
tially periodic potential. The interplay of the different 
time scales in the system, given by the period of the oscil- 
lating potential, the relaxation time of the deterministic 
dynamics and the diffusion time gives rise to non-trivial 
dynamical phenomena, such as an oscillation driven en- 
hancement of the effective diffusion. Otherwise if addi- 
tionally external forces are applied, the particle's motion 
becomes quasi-deterministic following the oscillations in 
direction of the applied force jumping several periods of 
the potential with minimal diffusion. Experimental sys- 
tems where such temporal modulations could be realized 
are free-flow diclcctrophoresis 18[ , colloidal particles in 
optical flclds [9 11], Josephson- junctions ll2Lll3l | or para- 



magnetic colloids in magnetic fields [IJ, |15[ . With small 
modification the dynamics can also describe neuronal ac- 
tivity being one type of a firing theta- neuron |16l. Il7|. 

We consider an overdamped Brownian particle under 
the infiucncc of thermal fiuctuations in spatially periodic 
potential with harmonically oscillating amplitudes. The 
corresponding Langcvin equation written in dimension- 
less form reads 



y = cos(y) sin(f2T) + F + V2D^{t), 



(1) 



As characteristic obscrvablcs we investigate the asymp- 
totic drift velocity and the effective diffusion coefficient 



Wdrift 
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with driving frequency fl, constant force F, and noise 
intensity D. 



Here (•) denotes the ensemble average. 

In the case of non- vanishing flux in the system (F ^ 0) 
we measure the quality of the directed transport by the 
so called Peclet number |18|: Pe = VdrHtL/Dcff where 
L is the characteristic length scale of the system. Here 
we set L equal to the spatial wave length A = 2tt. For 
Pe < 1 diffusion dominates the dynamics, and the di- 
rected transport plays a minor role in comparison with 
the non-directed spread of the probability distribution. 
For Pe > 1 the transport is dominated by the drift. The 
limit of Pe — >■ oo corresponds to a deterministic trans- 
port with vanishing effective diffusion on the character- 
istic length scale L. 

First we consider the case without bias in the system 
for F ~ 0. The spatial and temporal symmetry in the 
potential prevents any directed flux within the system 
("Adrift = 0). For large noise intensities the influence of the 
potential is negligible and the dynamics corresponds to a 
free Brownian motion with Z?off = D whereas in the limit 
of vanishing noise the interplay of two time scales controls 
the dynamics. On the one hand the external driving 
period T = 27r/i7 and secondly the intrinsic relaxation 
time Tr from an unstable potential maximum to a stable 
minimum. 

At oscillation frequencies much faster than the relax- 
ation time fl ^ 1/Tr the potential is self-averaged to be 
effectively fiat for the particle dynamics and the effective 
diffusion also converges to the free Brownian diffusion co- 
efficient D for n — )■ oo. However, for a relaxation time 
shorter than the half oscillation period t^ < T/2 = tt/H 
trajectories are able to approach the next minimum un- 
til the oscillation changes the minimum into a metastable 
maximum. Involving a small amount of noise the particle 
performs effectively discrete jumps between the minima 



of the potential for the extreme potential settings (max- 
imal barriers). 

This behavior can be easily described by an oscillation 
induced random walk. The equally distributed jumps to 
the left and right take place at discrete time intervals 
r/2 and correspond to the switching of the fixed points 
(two jumps per full period). The jump length is given 
by / = A/2 = TT. Thus the probability to find the par- 
ticle at the position yi = i ■ I after N jumps is given 
by the binomial distribution, which converges towards a 
Gaussian distribution in the long time limit describing 
diffusion with effective coefficient D^g = 7rO/2. As dis- 
cussed above for large frequencies the effective diffusion 
will reach asymptotically the noise intensity D. So the 
linear relation can not hold for the whole frequency range 
but has to pass a maximum. 

An initial Gaussian distribution around a stable fixed 
point (potential minimum) y„ splits as the fixed point be- 
comes unstable (transforms into a potential maximum). 
For intermediate frequencies only a certain fraction of the 
particle ensembles can reach the next minimum. The re- 
maining part moves back into the initial position. We 
attempt to describe this mechanism by approximating 
the probability u{^,D) of a particle starting at y„ + S 
{5 > 0) from which a particle is able to arrive at the 
neighboring fixed point ycnd = 2/n+i — £■ Here e repre- 
sents a small distance which the particle overcomes by 
fluctuations alone. This probability is given by the com- 
plementary error function 



u{n, D) = crfc 



(4) 



with the width u^ as the particle distribution at y„ and 
6 = 2arctan [exp (— 2/i7) cot (e/2)] as the initial devia- 
tion. The new expression for the effective diffusion coef- 
ficient has to be scaled by taking Eq. |4] into account and 
we get 
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The parameters e and cr^ remain unknown. As e repre- 
sents a minimal length and a^ is the variance of the initial 
Gaussian distribution, we make the following ansatz for 
their dependence on D: a^ ~ bD and e = e-\AD, where 
e, h are undetermined constants. This result is com- 
pared with numerical simulations. Examples are shown 
in Fig. [2] A reasonable choice of the two undetermined 
constants is e = 0.09 and b — f50. The linear increase 
of -Deff is recovered together with a maximum at ri^ax 
in agreement with the numerics (circles and squares in 
Fig. [2]). The inset in Fig. [2] shows the dependence of the 
position ilmax of the effective diffusion over the noise in- 
tensity D over three orders of magnitude. It illustrates 
that a single choice of e and b allows predictions on the 




FIG. 1: Effective diffusion vs. driving frequency from simu- 
lations for F = and D = 1.0 ■ 10"** (circles), D = 5.0 ■ 10"^ 
(squares) in comparison with results of of Eq. [5] (solid lines) 
for e = 0.09 and h — 150. The inset shows the prediction on 
the the position Slmax of the maximum of Deft vs. D (solid 
line) compared to numerical results. 



impact of noise on ilmax- Please note that the approx- 
imation is not valid at very large frequency values as 
limn-).oo -Deff (^, D) diverges. 

The addition of a constant force leads to a temporally 
constant tilt of the oscillating periodic potential. In this 
case a critical force can be defined, where for all times no 
potential barriers obstruct the drift motion of the parti- 
cle: i^crit = 1- In contrast to a static washboard potential 
we observe a finite drift for vanishing noise Z? — > also 
at subcritical forces F < Fcrit due to the oscillations of 
the potential. 

The drift speed Wdrift averaged over initial conditions 
2/ini S (0, 27r) as a function of the F shows different 
plateaus at f drift = 17, Sfi, 517, . . . , corresponding to 1:1, 
1:3 and 1:5 synchronizations (Figs. I2l3p . This behav- 
ior resembles characteristics of a driven oscillator, where 
such plateaus correspond to entrainment regimes of the 
oscillator to the external driving. Furthermore externally 
driven stochastic oscillators show a strong inhibition of 
the effective phase diffusion in the synchronized state [5[ , 
which is also observed in our system. 

Based on these similarities we attempt to describe the 
dynamics of our system in vicinity of the synchronization 
regime by a corresponding ansatz. We change into the 
CO- moving frame z{t) = y(t) — fit and consider the aver- 
aged deterministic dynamics over one oscillation period 
T = 27r/ri. Due to the undetermined time-dependence 
of z{t) we are not able to obtain a general solution of the 
averaged dynamics. We may write z as a Taylor series 
z{t') = z(0) + z'(0)t' + . . . with t' = t/T. Under the as- 
sumption z{t') changes slowly over one oscillation period 
we keep only the 0-th order and set z(t') = zt- With this 
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FIG. 2: Simulation of Eq. [T] at the 1:1 synchronization 
{Q — 5.0, D = 0.05): Deff (circles) and the Adrift (squares). 
The solid lines are the corresponding analytical results for the 
mean velocity (Stratonovich [I^]) and the effective diffusion 
from Eq. [3 



approximation wc obtain the stochastic Adlcr equation: 



(6) 



with A = \F-n\. 

Please note that due to the assumptions the equation 
holds only in the limit of fast oscillations with respect 
to the intrinsic relaxation time fl 3> t~^. Thus at 1:1 
synchronization [F = CI, Q ^ 1) we are in fact in the 
super-critical force regime. In the co-moving frame even 
for F > Fcrit — 1 the reduced Equation [6] describes a 
particle moving in a tilted stationary periodic potential. 

An analytical solution for Wdrift close to synchroniza- 
tion can be taken directly from the classical work of 
Stratonovich on the mean phase of an entrained stochas- 
tic oscillator |19[, whereas Doff can be calculated in the 
case of small noise D ^ 1/2, from the Kramers rates 
k to neighboring minima of the co-moving stationary 
potential I?cff = V (^^ + k~): 



Dcff — — cosh — n 
2tt \D 



,--J(Aarcsin(2A)+ii;o) 



(7) 



with Wo = iVl-4A2. 

The breakdown of effective diffusion in the synchro- 
nization regime correspond to giant Peclet numbers, 
which indicate a quasi-deterministic Brownian transport. 
For low noise strength {D < 0.01) and an optimal choice 
of F the system shows no dispersion even for extremely 
large times (i > lO^T, ~ 10'' particles). 

At low rj (< 1) the first synchronization regimes are 
located at subcritical forces F < 1. The dynamics can 
be described as a slip and stick motion dominated by 
the oscillating barriers. We assume a (S-peaked probabil- 
ity distribution at the begin of the sliding motion (slip 
phase) at the position tjq corresponding to the fixed point 
where the particle is located at the begin of the step. As 



soon as temporal barriers occur, the particle relaxes at 
yi = lX/2 with I = {2n+l) and n € N {stick phase). The 
envelope of the i5-like distributions of the relaxed particles 
is assumed to be Gaussian with a mean /Xg and width S. 
The probability to find the particle at the unstable fixed 
points y„i with m = 2n at the end of the step vanishes. 
The total probability to find the particle close to fixed 
point yi at the end of the step reads 



1 
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(8) 



Thus we can formally calculate the mean position and 
the variance after a single step 

CXD 

iVe) = ^ 2;2n-HrW2n-hl, (9) 

n— — oo 
oo 

{{y-{ye)f)^ Y^ {y2n+l~ {ye)fw2n+l, (10) 
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with y2n+i = {n + 5)A. As the distribution Pe{y) has 
a finite width and most of the probability will be con- 
centrated within few minima of the potential around 
the mean /ig. It is sufficient to consider only a fi- 
nite number of points n G [nmin , nmax] around fie 
and calculate the mean position and the variance by 
renormalizing the probabilities accordingly W2n+i = 
W2n+i/ X]n=n°^ u'2n-i-i- With rcsults obtained in Eqs. [9] 
and [To] we can calculate the mean drift and the effective 
diffusion as 



Adrift = 2{ye)/T, D,s = {{y - {ye)f)lT . 



(11) 



The mean of the particle probability distribution 
moves in the direction of the bias with an effective force 
G = G{F,n) within a time interval At(F, f2) caused 
by the combined effect of the constant bias F and the 
oscillating potential. The time of vanishing barriers is 
At — 2 aTCsm{F)/fl and for G we choose a linear depen- 
dence on F: G{F) = aF. Averaging over the potential 
oscillations leads to a = 1. 

The diffusion coefficient S in Eq. [8] depends on the 
applied force and shows a large increase close to the crit- 
ical force Fcrit hkcwise reported in [l| . For simplicity we 
assume here a linear increase S = bF. The simplified 
theory given in Eqs. [HI [11] is compared with simulations 
of the full system for fl = 0.03, 0.05 in Fig. [3] Sym- 
bols represent the numerical rcsults while the black solid 
line represents the analytical calculation for a = 0.8 and 
h = 0.04. The multiple peaks of the Peclet number with 
decreasing magnitude are reproduced together with the 
shift of the peaks towards larger F for increasing Q which 
leads to a decrease in the number of peaks for F < F^it ■ 
However the most obvious difference is the position and 
the height of the first peak of the Peclet number. 
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Via a transformation to the co-moving frame at fast os- 
cillations analytical expression based on previous results 
by Stratonovich and Kramers are found for the effective 
diffusion and mean velocity. At low frequencies we de- 
scribe the system dynamics through a simplified ansatz 
of the transport process, which for reasonable choice of 
effective parameters are in agreement with the numerical 
results for < i^ < Fcrit- Our results suggest that in 
the unbiased case the introduction of an spatio-temporal 
periodic force enhance the diffusive transport, whereas 
for a finite bias it may significantly improve the quality 
of directed Brownian transport by suppressing thermal 
fluctuations. 
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FIG. 3: Comparison of simulations (symbols) with the the- 
oretical results at low fi (solid lines): mean velocity Udrift, 
effective diffusion Deff and Peclet number Pe vs. bias F for 
SI = 0.05 (left column) and SI = 0.1 (right column). The 
black lines show the linear ansatz G = aF with a = 0.8, the 
red lines show the non-linear ansatz with G = a{F)F with 
ai = 1.6, a2 = 0.8, c = 10 for both frequencies. The only 
changing parameter is Ft: for Q = 0.05 Ft — 0.28, for fi = 0.1 
Ft = 0.4. For both examples shown b = 0.08. 



A small bias F is potentiated by the slope of 
the oscillating potential therefore the parameter a 
should be chosen as a > 1. We describe this 
through a step function around Ft as '-aiF) = 02 -t- 
(01 — 02) [tanh {c{Ft — F)) — 1] /2 This phenomenological 
modification of the model leads to the red (gray) line and 
shows agreement to the numerical simulations. Please 
note that Ft is the only changing parameter in calcula- 
tions of the non-linear ansatz for different frequencies. 
This good agreement of the theoretical results with the 
numerical simulations shows that the theory accounts for 
the decisive mechanism responsible for the observed be- 
havior: the combination of transport (slip phase) with re- 
peated confinement (stick phase) at temporal minimum 
of the potential. 

In this work we have analyzed the impact of a spatio- 
temporal oscillating potential on Brownian transport. 
We have demonstrated how the oscillating potential may 
enhance or suppress the effective diffusion. We derived 
an expression for Des (fi, D) for vanishing bias F = 0. We 
have obtained an optimal driving frequency which max- 
imizes the diffusive transport. For finite bias _F > we 
have shown the occurrence of synchronization regimes, 
with strongly suppressed effective diffusion and quasi- 
deterministic transport at finite noise strengths revealed 
through giant Peclet numbers. 
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